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p • Abstract 

D ; 

^fn ■ We describe the calculation of inclusive Higgs boson production at hadronic 

colliders at next-to-next-to-leading order (NNLO) in perturbative quantum 
^ ! chromodynamics. We have used the technique developed in reference [4]. 

H ' Our results agree with those published earlier in the literature. 



1 Introduction 

The discovery of the Higgs boson will shed light on the symmetry break- 
ing machanism of the Standard Model (SM). The experimental bound from 
the LEP experiments and precision studies within and beyond the standard 
model strongly suggest that hadron colliders such as the Tevatron and the 
LHC will see the Higgs boson if it exists. At these machines the dominant 
contribution to single Higgs boson production is the gluon-gluon fusion pro- 
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cess through heavy quark loops. The reason for this is that the Higgs boson 
couples strongly to heavy quarks. In addition the large gluon flux at the 
LHC enhances the total inclusive cross section substantially. The NLO cor- 
rections along with NLO parton distributions yield a large K factor and also 
show a strong scale dependence. Hence there is need for improved parton 
densities using NNLO splitting functions as well as inclusion of NNLO par- 
tonic cross sections. The NNLO correction to Higgs boson production was 
first computed by an expansion technique in [1]. Exact results were obtained 
in [2] using Cutkosky rules. In our work [3], we have used a straightfor- 
ward technique which was adopted in [4] to compute the NNLO corrections 
to Drell-Yan process. Clever choice of the integration variables in specific 
frames makes the computation manageable. 

2 Method of Computation 

We use the effective Lagrangian approach which emerges from the SM in the 
heavy top quark limit {rrit oo) and is found to be a good approximation 
at hadron colliders. For Higgs boson production at hadron colliders at the 
NNLO level one has to compute 1) tree level a + b c + d + H, 2) one-loop 
corrected a + b ^ c + H and 3) two-loop corrected g + g ^ H, where a, b, c, d 
are light partons such as quarks, antiquarks and gluons whose interactions are 
governed by QCD. These corrections involve the computation of 2 — 3 body 
phase-space integrals and two and one-loop momentum integrals followed by 
2 — > 1 and 2^2 phase-space integrations respectively. We use dimensional 
regularization (space time dimension is taken to be n = A+e) to regulate both 
ultraviolet and infrared (soft and coUinear) divergences. We first describe 
here how we have performed three-body phase-space integrals for 2^3 tree 
level matrix elements and the two-body phase-space integrals for one-loop 
corrected matrix elements. The 2^3 body processes involve two angular 
integrations (say 9,(j)) and two parametric integrations (say z,y). Before we 
perform these integrations, it is important to classify the matrix elements 
in such a way that the phase-space integrations over them can be done in 
suitable frames. For example, when the Higgs boson is produced from the 
incoming partons, the center of mass (CM) frame of outgoing partons is the 
most suitable frame, because the massive propagators 1/(P{|P25) {c(,P ^ 1) 
will not involve angular dependence in this frame. Here Pi5 — (pi + puf 
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where ph is the momentum of the Higgs boson, and pi is the momentum 
of the massless parton. Similarly when the Higgs boson is produced from 
an outgoing parton, we choose the CM frame of incoming partons where 
1/(P^P^) (cK, /3 > 1) do not depend on the angles. Complications arise when 
we encounter processes where the Higgs boson is produced by both initial 
and final state partons, i.e., where interference terms of the form 1/(P{|P^^) 
{a,j3 > 0) appear. In this case we have chosen the CM frame of the 4th 
(or 3rd) parton and the Higgs boson where the angular integrals and other 
parametric integrals are less difficult. In the CM frame of the incoming 
partons and the CM frame of the outgoing partons we perform the angular 
integrations exactly using the result given in [5]. 



where = (1 — cos^)~*(l — cos x cos ^ — sin % cos sin ^)~-^. Here cos^ is 
related to kinematical variables such as x (= m'j^/s, s-CM energy) and the 
integration variables z and y, and Fi^2 is the hypergeometric function. In 
the CM frame of the 4th parton and Higgs boson, due to the complexity, 
we performed one angular integration (say 0) exactly and performed the 
remaining 6 integration after expanding the integrands in powers of £ = n— 4. 
In all these frames, using various Kummer's relations, the hypergeometric 
functions are simplified to the form Fi^2(±£/2, ±£/2, lie/2; T>{x, y, z)) which 
is the most suitable for integrations over z and y. The next hurdle in the 
computation is the appearance of terms with large powers in 1/(1 — z) or 
1/z. We have reduced the higher powers of 1/(1 — ^)"+^^ or 1/2;"+^^ where 
a > 1 by successive integration by parts with exact hypergeometric functions 
until we arrive at 1/(1 — zy~^^^ or l/z^^^"^ multiplied by regular functions. 
We have used the following identity to accomplish this: 
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r de r d(f) sin"-^ e sin"-^ (f) Cij {e, 0, x) = 2^~'~^' 
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l;cos2(x/2)). 



dz 




(2.1) 
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In the end we are left with integrations of the form Jq dzz ^ ^^f{z) and/or 
Jq dz{l — z)~^~^^ f{z). Such integrals are simplified as follows: 

dzz-^-^'f{z) = dz z-^-^'f{z) + dz z-^-^'f{z) 5«l. (2.2) 
Jo Jo Js 

The first term can be evaluated to be f{Q){[l3eY^ + \og5+[pe/2] log^ 5 + ■ ■ ■). 
After expanding z~^^ in powers of e in the second term the z integration can 
be performed exactly order-by-order in e with non-zero 5. At the end the 
5 dependence cancels in each order in e. Since the z integration over the 
hypergeometric functions is nontrivial due to their complicated arguments, 
we have expanded them in powers of e prior to the z integration: 

Fi,,(|,|,l + |;P(^)) = l + ^Li2(P(^)) + ^(Si,2(I^(^))-Li3(l^(^))) 

+ |^(Si,3(^^(^)) - ^2,2{V{z)) + U,{V{z))) (2.3) 

where S„,p(^) = (-l)"+P-i[(n - l)!p!]-i /o c^^M"^ log"-^(t) log^(l - zt) with 
n,p > 1 and Lin{z) = S„_i,i(2;)(see [6]). We have repeated the same proce- 
dure to perform the remaining y integration. In addition to 2 3 contri- 
butions, we encounter one-loop corrected 2^2 processes at NNLO level. 
Here the one-loop tensorial integrals are reduced to scalar integrals using the 
Passarino-Veltman reduction procedure implemented in n dimensions. The 
resulting one-loop two and three-point scalar integrals can be expressed in 
terms of kinematic invariants. In the case of the four-point function, the 
scalar integrals can be expressed only in terms of hypergeometric functions 
which increases the complexity of the two-body phase-space integrations. 
We follow the procedure adopted for the 2 — >■ 3 phase-space integrations to 
perform the two-body phase-space integrations. After performing all these 
integrals, we have removed all the ultraviolet divergences by strong coupling 
and operator renormalization constants. The remaining collinear divergences 
are removed by mass factorization. Then we are left with finite partonic 
cross sections which are folded with parton distribution functions to com- 
pute hadronic cross section for the inclusive Higgs boson production. 

Using the method described above we have successfully computed the 
NNLO corrections to Higgs boson production at hadron colliders and found 
complete agreement with the results of [1, 2]. We find that NNLO corrections 
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improve the convergence of the perturbative result and decrease the scale 
ambiguities inherent in it. 
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